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We propose a scheme for entangling the motional mode of a trapped atom with a propagating 
light field via a cavity-mediated parametric interaction. We then show that if this light field is 
subsequently coupled to a second distant atom via a cavity-mediated linear-mixing interaction, it is 
possible to transfer the entanglement from the light beam to the motional mode of the second atom 
to create an EPR-type entangled state of the positions and momenta of two distantly-separated 
atoms. 
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PACS numbers: 03.67.Hk, 32.80.Lg, 42.50.-p 



I. INTRODUCTION 

The generation, distribution, and application of con- 
tinuous variable quantum entanglement are topics of con- 
siderable interest at present, spurred on in large part 
by the burgeoning fields of quantum communication and 
quantum computation. In this context, a variety of pro- 
tocols for continuous quantum variables have been pro- 
posed and in some cases already demonstrated, including 
uantum teleportation ||, ||, quantum cryptography 
Hi H 01 > quantum dense coding [g[ ^, ^ , and quantum 
computation 

These protocols have in large part focussed on im- 
plementations involving nonlinear optics and propagat- 
ing light fields, which has obvious advantages in terms 
of long-distance communication and existing quantum- 
optical technology. However, for purposes of storage (i.e., 
memory) and local manipulation, a number of alternative 
physical systems are being actively investigated and ap- 
pear very promising; notabh;, collective atomic spin sys- 
tems |l5], |l6, 17, ISt and quantized vibrational 
states of trapped atoms ]19|, |2C|, pT| , p2| , . Of particular 
interest is the capability of establishing long-lived entan- 
glement of the Einstein-Podolsky- Rosen (EPR) type ||2j] 
between actual or effective position and momentum vari- 
ables of two separated atomic systems. This capability 
has in fact been demonstrated recently in an experiment 
involving collective spins of a pair of atomic ensembles 
and nonlocal Bell measurements using off-resonant light 
pulses 0. 

A scheme for preparing an EPR-type state of the ac- 
tual positions and momenta of a pair of atoms has been 
put forward in [^ , making use of interactions in cavity 
quantum electrodynamics (cavity QED) to facilitate the 
transfer of quantum correlations from light fields to mo- 
tional modes of tightly trapped atoms. The source of the 
quantum-correlated light fields was taken to be an optical 
nondegenerate parametric amplifier (see, e.g., j2^, [26|). 

Here, we present an alternative approach to preparing 
such a motional state which, in contrast to the scheme 
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of IgOf, does not require a separate source of quantum- 
correlated light beams. Furthermore, unlike a number 
of other proposals, it does not require entangling mea- 
surements to be made, or a carefully timed sequence of 
suitably shaped light pulses. Through an atom-cavity 
coupling similar to, but modified from that of |20| , an 
effective parametric interaction between cavity and mo- 
tional modes generates continuous variable entanglement 
between the motion of the trapped atom and the light 
field exiting the cavity. The entanglement "carried" by 
this light field can subsequently be distributed to a dis- 
tant location and thence to another atom (or atoms). 



II. THE MODEL 

The essential details of the scheme to be described 
in this work are illustrated in Fig. 1. A pair of atoms 
are harmonically confined inside separate optical cavities, 
with the light exiting one of the cavities coupled into the 
second cavity (but not vice- versa). Auxiliary lasers, inci- 
dent through the sides of the cavities, combine with the 
cavity fields to drive Raman transitions between neigh- 
bouring vibrational levels of the motion of each atom. 
Note that only a single internal atomic state (i.e., a sta- 
ble ground electronic state) is assumed to be involved, as 
will be discussed below. 



A. Motion-Light Coupling 

The laser-atom-cavity interactions responsible for the 
coupling between motional and light modes have been 
discussed in detail previously (see, e.g., |2^), but for 
completeness we include a brief description. Considering 
just a single-atom configuration, our system is modeled, 
in a frame rotating at the laser frequency wl, by the 
Hamiltonian 

j=x,y,z 

+ TT'[£i.{y,z,t)a+ + £l{y,z,t)aJ\ 
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FIG. 1: Schematic diagram of the system. Atoms inside the 
cavities are tightly confined in harmonic traps. The cavity 
and laser frequencies are chosen to drive Raman transitions 
between neighbouring vibrational levels of the motion. 
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The first line describes the quantized harmonic motion of 
the trapped atom, with bj the annihilation operator and 
i/j the vibrational frequency for motion along the j-axis. 
The quantized cavity mode, with annihilation operator a 
and frequency Wcav, is detuned from the laser frequency 
hy S = LUcnv — i-^L- The ground and excited electronic 
states of the atom that are coupled by the light fields, \g) 
and |e), are separated in energy by hujo, and the detuning 
of the atomic transition from the laser frequency is given 
by A = LJo — ^L', tT+ — \e){g\ and cr_ = \g){e\ are the 
atomic raising and lowering operators, respectively. 

The laser field is treated as a classical field of (complex) 
amplitude £h{y, z, t) and is assumed to propagate in the 
yz-plane. The cavity mode is aligned along the x-axis 
and its coupling to the atomic transition is described by 
the last line in (yj) , where go is the single-photon coupling 
strength and k is the wavenumber of the cavity field. The 
choice of a sine function, with x = [T%/2mVx)'^^'^{bx + ^1) 
the position operator of the atom, denotes that the trap is 
assumed to be centered at a node of the cavity standing- 
wave field. 

By taking the laser-atom detuning A to be large (i.e., 
A ^ {vx, 5, ffOj 7}: where 7 is the spontaneous emis- 
sion linewidth of the state |e)), population of the excited 
internal state |e), and hence spontaneous emission, can 
be made negligible. From the Heisenberg equation of 
motion for the atomic lowering operator, 

(T_ = — iAcT- + iShiu, z, t)az + igo shi{kx)aaz, (2) 

one can then take 

-£i.{y,z,t) gosm{kx)a 
- X X ' (3) 



noting that Uz — and CFz\g) = —\g) (i.e., setting CTz ~ 
— 1 in the equation of motion for u- ) . 

Now, if the size of the harmonic trap is small compared 
to the optical wavelength (Lamb-Dicke regime), then we 
can make the approximation 



sin(fca;) ~ r-jxib^ + 6^), 



(4) 



where 77^; <C 1 is the Lamb-Dicke parameter. In this case, 
one can also assume a configuration such that the posi- 
tion dependence of the laser amplitude over the extent of 
the trap can be ignored, i.e.. 



£h{v,z,t) ~ £i,{t)e 



(5) 



and hence only motion along the x-axis is of relevance. 

Finally, for the regimes we are most interested in, the 
cavity mode is only ever weakly excited. This fact, com- 
bined with the smallness of the Lamb-Dicke parameter, 
enables us to neglect the second term in (^ in compari- 
son to the first term, reducing our approximate solution 
for (T_ to the simple form 



fL(t)e 



(6) 



Using this form, and shifting the zero of energy to remove 
scalar terms, the Hamiltonian describing the cavity mode 
and motion along the a;-axis reduces to [pTl 
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Under suitable conditions and with appropriate choices 
of the cavity-laser detuning, b = Wcav — t^L, we are able to 
choose between a parametric or linear mixing interaction 
between the cavity and motional modes, as we now show. 



B. Cavity 1: Parametric Amplification 

For convenience, we now drop the subscripts x and L 
and use the subscript 1(2) to denote system 1(2). For 
cavity 1 we choose <5i — —vi- Moving to an interaction 
picture to remove the systematic motion associated with 
the first two terms of (^, we obtain 



Hi = ni{t) [aibie"^' + a\ble 



where 



- a\bie 



2iy-it—i(f>i 



■ aib{e 



t „2ifit+i0i 



hrjigiSi (t) 
Ai ■ 



(8) 



(9) 



If we now assume that the trap frequency i^i is large, 
such that i>i 3> |r2i(t)|, and also i>i 3> Ki, where ki is the 
amplitude decay rate of the cavity field, then the rapidly 
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oscillating terms in (^) can be neglected in a rotating- 
wave approximation to give 



(10) 



This is the Hamiltonian for a nondegenerate paramet- 
ric amplification process and via this process continuous 
variable entanglement can be generated between the mo- 
tional and cavity modes in system 1. 



C. Cavity 2: Linear Mixing 

For cavity 2 we choose 82 = 1^2- Moving to the ap- 
propriate interaction picture and employing the rotating- 
wave approximation once again under the condition that 
1^2 ^ 1^2(01 ^2 ^ K2, we derive for system 2 the 
effective Hamiltonian 

H2 = n2{t) (4626^*^^ + bla2e"^^) , (11) 



where 



n2it) 



(12) 



This describes a linear mixing interaction which, as we 
shall see, can lead to an exchange of properties between 
the cavity and motional modes. 



D. Cascaded System 

Having specified the effective interactions occurring in 
each atom-cavity system, it should now be quite clear 
what our aim is. Via the parametric interaction, light in 
cavity mode 1 is entangled with motional mode 1. This 
light ultimately exits cavity 1 due to cavity decay at the 
rate ki and it may then be coupled into cavity 2, where, 
through the linear mixing process, entanglement of the 
light field with motional mode 1 may be transferred to 
motional mode 2, thereby entangling the two motional 
modes. 

To examine this transfer process in detail, we turn nat- 
urally to a cascaded systems model |9[ p2| , 
which enables us to describe a unidirectional driving of 
cavity 2 with the output light from cavity 1. If we assume 
that the dominant input and output channel to and from 
each cavity is through just one of the two mirrors form- 
ing each cavity (as depicted in Fig. 1), then the master 
equation for our cascaded system is 

+ Ki{2aipa\ — a\aip — pa\ai) 
+ K2{2a2pa\~ a\a2P ~ pa\a2) 
- 2^/eK^{[a\,aip] + [pa\,a2]), (13) 



where 

H = Hi+H2, (14) 

and the parameter e, satisfying < e < 1, accounts for 
losses in transmission and for coupling inefficiency. Ideal 
transmission and coupling corresponds to e = 1. Note 
that implicit in the form (13) is the assumption that the 
two cavity frequencies are equal, and hence that the dif- 
ference in laser frequencies is set to 

^Ll - ^L2 = i^l + l'2- (15) 



III. ADIABATIC APPROXIMATION 

In the (overdamped) regime where both ki and K2 are 
much larger than the coupling rates f2i and ^2^ our model 
can be simplified further by adiabatically eliminating the 
cavity modes from the dynamics. In particular, the mas- 
ter equation (4) can be written in the form 



p = (£0 + l^c)p, 



where 



— 1 

Cap = —[H,p\, 



(16) 



(17) 



CcP = Hii{2aipa\ — a\aip ~ pa\ai) 

+ K2{2a2pa\- a\a2P - pa\a2) 

- 2^eKiK2{[a\,aip\ -|- [pa\,a2])- (18) 

In the adiabatic limit, a master equation for the reduced 
density operator, pb, of the motional modes alone can be 
derived as (see, e.g., [p2[) 



Ph = Trc <^ £1 



dr f 



i^OPl \ Pb 



(19) 



where the trace is taken over the cavity modes and 
is the steady state density operator of the cavity modes 
satisfying CcP% — 0. 

Evaluation of (^9|) requires calculation of steady-state 
two-time correlation functions of the cavity operators. 
From the master equation pc — CcPc one obtains the 
following solutions for the mean cavity amplitudes (for 
t > 0) 

(ai(t))c = e-«i*(ai(0)) 
(a2(t))c - e-''i*(a2(0)) 

2^KlK2e 



(20) 



-K2t 



'^^*)(ai(0)),.(21) 



K2 — Kl 

With {aia\)l — {a2al)l = 1 being the only nonzero 
steady-state equal-time correlations, the quantum regres- 
sion theorem gives 

{ai{T)a\{0))l = {ai{0)a\{T))l = e-^^\ (22) 

(a2(r)4(0))? = (a2(0)4M>c = e-«^^ (23) 

(a2(r)4(0))^ ^ (ai(0)4M>c 

_ 2^KiK2e 



K2 — Kl 



1, (24) 
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as the only non-zero two-time correlation functions (re- 
sulting from the master equation = CcPc)- 

Using these correlation functions in (|l^), the master 
equation for the reduced density operator of the motional 
modes is 



Pb 



Ti{2b\pbhi - bib\pb ~ Pbbib\) 
^2{'2b2Pbbl - 4^2/56 - Pbb\b2) 



2^eTiT2([blb\ 



Pb\e 



-i(0i-02) 



where 



[Pbbi,b2]e 



Ti = ^ and T2 = — 



i(0l-02) 



K2 



(25) 



(26) 



are the effective growth and decay rates, respectively, for 
motional modes 1 and 2. In this reduced model a direct 
coupling between these modes now appears in the form 
of the last term in (ESf). 



IV. MOTIONAL STATE ENTANGLEMENT 

A. Motional Mode Correlations 

From the master equation (|25|), a closed set of (inho- 
mogeneous) differential equations can be obtained for the 
correlation functions {b\bi), (62^2) and (6162)- These are 
(setting 01 = 02 for simplicity) 



(61&2) = (ri-r2)(foi62) 

+2 VeriFa {{b\b,) 
{b\b,) = 2T,{l + {b\b,)), 
{bib2) = -2T2{blb2) 



(27) 
(28) 



-2^er^{{b\bl) + (6162)), (29) 



with solutions 

{b\mit)) 

{bi{t)b2{t)) 

(bimit)) 



2\/eTlT2 Tit/ Fit 



ri + r2 
4£rir2 
(ri + r2 



-r2t\ 



r2t\2 



(30) 
(31) 

(32) 



where both oscillators are assumed to have initially been 
in their ground states. Note that, under these conditions, 
(6i(<)6i(i)) = {b2{t)b2{t)) = {b\{t)b2{t)) = for all t. 

The growing value of (6i(i)&2(^)) demonstrates that 
our cascaded system enables correlations to develop be- 
tween the two motional modes. The nature of the cor- 
relations is, not surprisingly, reminiscent of a two-mode 



squeezed state and leads us to examine correlations be- 
tween the positions and momenta of the trapped atoms. 

Before proceeding, however, we note that the exponen- 
tial growth exhibited by the mean excitation numbers 
must eventually become problematical for our model, 
since the Lamb-Dicke assumption breaks down once the 
excitation numbers become large enough that the phys- 
ical extent of the atomic wavepacket is no longer much 
smaller than the wavelength of the light. We shall return 
to this point when we consider practical aspects of the 
scheme. 



B. Position and Momentum Variances 

Defining position and momentum operators for the mo- 
tional modes as 

Xj^b,+b], P,=^i{b,-b]), (j = l,2), (33) 

variances in the sum and difference operators are calcu- 
lated to be 



{{x,±x2f)^mTP2f 



^rit , VfTir^ Pit 

^ ri + r2^' 



-,-^2t^ 



(34) 



where the "vacuum" level (i.e., where the atoms are both 
in their ground motional states) is 2. We note first that, 
since (e^i'-e-r^*) > 0, {{X^+X2?) = {{P1-P2?) > 2 
at all times, i.e., these variances are only ever increased 
compared to the vacuum level. 

However, the variances ({Xi — X2Y) and ((Pi + ^2)^) 
can be reduced below the value of 2, as illustrated in 
Fig. 2, where {{Xi— X2Y) is plotted as a function of time 
for various values of the ratio r2/ri, and with e = 1. The 
case r2/ri = 1, where the effective damping rates of the 
two motional modes are equal, is particularly interesting 
and leads to the simplified result 



((^1 - X2 



= {{Pl+P2f) 

^ 2e"2^i* ^ 



as t 00. 



(35) 



So, using this scheme it is, in principle, possible to pro- 
duce a state corresponding to perfectly correlated posi- 
tions and perfectly anti-correlated momenta of the two 
atoms, i.e., to produce an EPR state. Once again, in 
comparing this approach to that of it should be 
emphasized that the present scheme does not require a 
separate source of entangled light beams. 

If r2/ri 7^ 1 and/or e < 1, then the variance {{Xi — 
^2)^^) attains a finite minimum value at a particular time, 
after which it increases indefinitely. Defining A = r2/ri, 
this minimum value is calculated to be 



FIG. 2: The variance {{Xi — X2)^) as a function of time for 
Fi/ri = 0.2 (pluses), r2/ri = 0.5 (dot-dashedline), r2/ri = 
1 (solid lino), r2/ri = 2 (dotted line), and r2/Ti = 3 (dashed 
line). The coupling is taken to be ideal, i.e., e = 1. 



{{Xi — X2)^)min — 



{1 + X - 2^/Xi)^h 



-A 
1 + A 



l + A-2\/Ae, 



(36) 



and the time at which the minimum value occurs is given 

by 



rifmin(A) 



1 



1 + A 



■In 



2\VXe 



l + X-2y/X^ 



(37) 



Plots of ((Xi — X2)^)niin and Fi^min vBrsus A are shown 
in Figs. 3 and 4 for several values of e. Note that imin > 
only where A > l/(4e). If A < l/(4e) then there is no 
reduction in the variance below the vacuum level. 

From Figs. 3 and 4 (in particular, from their asymme- 
try about the point A = 1), it is clear that the scheme 
generally performs best in the regime where r2 > Fi (i.e., 
A > 1). In particular, significant reductions in the vari- 
ances below the vacuum level occur over a broad region 
of r2 values, provided r2 > Fi. 

With decreasing values of the parameter e (i.e., with 
decreasing coupling efficiency or increasing transmission 
losses), the minimum attainable variance increases and 
occurs for larger values of A and at somewhat earlier 
times. It follows that, for a realistic system with e < 1, it 
would be advantageous to work in a regime with r2 > Fi. 
This is further illustrated in Fig. 5, where the variance is 
again plotted as a function of time, but now for several 
combinations of A and e. It is worth noting that a signif- 
icant reduction in the variance (^ 73%) is still possible 
with e = 0.8, i.e., ((Xi - X2)^)i„in 0.54 at Fifmin ^ 0.8 
for A = 2. 




FIG. 3: Minimum variance {{Xi — X2)'^)min as a function 
of A = r2/ri for e = 1 (solid line), e = 0.9 (dashed line) 
(minimum value of 0.30 at A = 1.46), and e = 0.8 (dot-dashed 
line) (minimum value of 0.54 at A = 1.95). 



It should also be noted that, for the examples chosen in 

the regime where e < 1 and A > 1 , the minimum variance 
generally occurs at times Fit of the order of 1. The mean 
excitation number for motional mode 1, {b\bi), is then of 
the order of — 1 ~ 6, with a slightly smaller value for 
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FIG. 4: Time at which minimum variance occurs, Titmin, as 
a function of A = r2/Ti for e — 1 (sohd hne), e = 0.9 (dashed 
line), and e — 0.8 (dot-dashed hne). 



reduced and the dynamics obviously become somewhat 
more complicated, e.g., oscillatory behaviour starts to 
feature and excitation of the cavity modes increases. 
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FIG. 6: Comparison of numerical and adiabatic approxima- 
tion results for {{Xi — X2Y) versus time with ki = K2 = 1. 
The adiabatic approximation is shown in each plot as a 
solid line, (a) A = 1, with f2i = ^2 = 0.1 (dashed line), 
Hi = 02 = 0.2 (dot-dashed hne), ^1=^2= 0.5 (dotted 
line), (b) A = 2, with fii = 0.1, 0,2 = \/2 0.1 (dashed 
line), 0.1 = 0.2, O2 = ^2 x 0.2 (dot-dashed line), fii = 0.5, 
0.2^ V^x 0.5 (dotted line). 




FIG. 5: Variance ((Xi - Xa)^) as a function of time for (a) 
A = 1 and (b) A = 2, with e = 1 (solid line), e = 0.9 (dashed 
line), and e = 0.8 (dot-dashed line). 



V. COMPARISON WITH NUMERICAL 
CALCULATIONS 

Starting from the master equation ( [l3| ) , which includes 
the cavity dynamics, a closed set of differential equations 
can be obtained for various correlation functions of the 
system. We have integrated these equations numerically 
(using a 4th-order Runge-Kutta method) and computed 
the variance ((^1 — X2f) for comparison with the re- 
sults of the adiabatic approximation. This comparison 
is presented in Figs. 6 and 7 for several example sets of 
parameters. 

On the timescales shown, the adiabatic approximation 
is seen to work well provided the coupling strengths 
are, roughly, at least five to ten times smaller than the 
cavity decay rates Kj. For larger values of fij the maxi- 
mum degree of reduction in the variance is significantly 
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FIG. 7: Comparison of numerical and adiabatic approxima- 
tion results for {{X\ — X2Y') versus time with ki = ^2 = 1 
and nonideal coupling between the cavities (i.e., e < 1). The 
adiabatic approximation is shown in each plot as a solid line, 
(a) A = 1, with r2i = ^2 = 0.1 (dashed lines), (b) A = 2, with 
fii = 0.1, ^2 = ^2 X 0.1 (dashed lines). 



VI. PRACTICAL CONSIDERATIONS 

We now consider in slightly more detail the conditions 
under which the most significant assumptions required by 
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our model should be satisfied. Firstly though, as a more 
general comment, we note that exciting progress has been 
made recently in experimental cavity QED with single 
trapped atoms or ions In fact, various 

"ingredients" of the scheme presented in this work have 
already been demonstrated. 



A. Trap Frequency 

The neglect of terms in the effective motion-cavity 
mode interaction Hamiltonians which vary like e^^*'^'* 
requires that the trap frequencies be large in comparison 
with the cavity decay rates Kj and the effective coupling 
parameters flj. To quantify this a little more precisely, 
previous numerical studies have shown that this rotating- 
wave approximation is very good provided the trap fre- 
quencies are at least an order of magnitude larger than 
Kj and Qj |2^. We note that an experimental situ- 
ation with h'j 3> K has been realized recently with single 
calcium ions trapped inside a high-finesse optical cavity 



B. Lamb-Dicke Approximation 

The Lamb-Dicke approximation has been examined in 
some detail in |^ . Taking into account the spread of the 
phonon number distribution associated with a general 
state, a condition for the validity of this approximation 
can be derived as 



1 



(38) 



where is the mean phonon number, cr? is the variance 
of the number state distribution, and a ~ 2 — 3 (i.e., a 
few standard deviations from the mean). If we assume 
that the number state distribution of each mode is close 
to that of a thermal mode, then we can take — (n^ + 
fixY^'^ ~ fia; -|- 1/2 for fix > 2, and, with a = 3, the 
condition becomes 



1 
2^ 



(39) 



For rjx = 0.1 this reduces to ^ 49. As noted ear- 
lier, at the time when minimum variances are generally 
achieved by our scheme (e.g., with e < 1 and A > 1), the 
mean phonon excitation numbers are of the order of 5-6, 
suggesting that Lamb-Dicke parameters of the order of 
0.1 or smaller are sufficient. 

Such values of the Lamb-Dicke parameter have been 
achieved with single atoms in cavity QED settings ||3^, |3^, 
|37| . An ability to position the trap center very precisely 
at any point along the cavity standing wave field (e.g., 
at a node of the field, as required by the present scheme) 
has also been demonstrated p7[ . 



Before continuing, we briefly mention the further in- 
teresting possibility of introducing a nontrivial time de- 
pendence to one or both of the coupling laser fields. One 
such example is illustrated in Fig. 8. Here, SI2 = 1 (in 
dimensionless units) is fixed, while fli{t) — fl2Sm{t/T) 
(r — 20,25,30), i.e., the effective parametric coupling 
strength is slowly increased from zero in the first atom- 
cavity system. If we consider the point for each case at 
which the variance is reduced to 0.2 (i.e., a 90% reduction 
below the "vacuum" level), then the corresponding mean 
phonon number at this time is seen to be only of the 
order of 2-3 for time-dependent ili(t), whereas for the 
case of constant Jli and ^2 (also shown in the figures) 
a variance of 0.2 is only attained with a mean phonon 
number greater than 10. 

We have not explored other time dependencies in de- 
tail, but it would appear that there could be some advan- 
tage to doing so, in particular, from the point of view of 
minimizing the variance while maintaining mean phonon 
excitation numbers which comfortably satisfy the Lamb- 
Dicke constraint. 
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FIG. 8: Variance {{Xi - Xif) and mean phonon number 
{b\bi) as a function of time from numerical simulation of 
Eq. ^ using e = 1, ki = K2 = 10, ^2 = 1 (r2 = 
n|/K2 = 0.1), and Qi{t) = ^2 sin(t/20) (dot-dashed line), 
Q.i{t) = ^2 sin(t/25) (dashed line), Sli(f) = ^2 sin(t/30) (dot- 
ted Une). The case £7i = n2 = 1 is shown as a solid line. 



C. Atomic Spontaneous Emission 

In practice, finite excitation of the atomic excited state 
will introduce the effects of spontaneous emission to the 
dynamics. Random momentum recoils associated with 
spontaneous emission events will disturb the motional 
states in an uncontrollable manner. However, if the rate 
at which such events occur is small compared with the 
rate at which the entangled states of interest are prepared 
by our scheme, then the effects of spontaneous emission 
can essentially be neglected. 
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In the rate of spontaneous emission events affect- 
ing the motional state is estimated for the kind of con- 
figuration we are considering here. By comparing this 
rate with the rates F^, which characterize the prepara- 
tion of the entangled motional states, the general condi- 
tion under which the effects of spontaneous emission can 
be neglected is derived to be 



» 1, 



(40) 



where 7 is the linewidth (FWHM) of the internal atomic 
excited state. This is basically the condition of strong- 
coupling cavity QED, as realized, for example, in the 
single-atom experiments of |35| (where values of 

go/{Kj) in the range 30-150 were achieved). 

Finally, note that throughout this work we have ne- 
glected any forms of motional state decoherence asso- 
ciated with the trap itself. This seems reasonable, at 
least in the case of trapped ions, where typical timescales 
for motional decoherence and heating observed in recent 
experiments are of the order of milliseconds or longer 
p8| , 3^. With values of go and k in the MHz range, we 
would anticipate a timescale for preparation of the mo- 
tional states (i.e., F^^^) that is more than an order of 
magnitude shorter than such decoherence times. 



metric (or two-mode squeezing) interaction realized in a 
single-atom cavity QED setting. Through decay of the 
cavity field through one of its mirrors, the cavity field 
entanglement is transferred to an external propagating 
field which may be coupled to the motional mode of a 
distant atom via a second cavity-QED-mediated inter- 
action. This enables the preparation of an EPR-type 
entangled state of the position and momentum variables 
of separated atoms. This has a distinct advantage over a 
related scheme for preparing such a state in that a 
separate source of entangled light beams is not required. 

We have given consideration to various practical is- 
sues associated with the scheme and pointed to a variety 
of recent experiments in single-atom cavity QED which 
collectively offer great encouragement to our proposal. 
Realization of the scheme would offer the exciting pos- 
sibility of implementing a variety of continuous variable 
quantum computation and communication protocols, as 
well as opening the door to further investigations of fun- 
damental aspects of entangled quantum systems, only 
now with distantly-separated massive particles [EOl . 
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